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ABSTRACT 

We analyze the uncertainties and possible systematics associated with the 
"Dark Flow" measurements using the cumulative Sunyaev-Zeldovich effect com- 
bined with all-sky catalogs of clusters of galaxies. Filtering of all-sky cosmic 
microwave background (CMB) maps is required to remove the intrinsic cosmo- 
logical signal down to the limit imposed by cosmic variance. Contributions to 
the errors come from the remaining cosmological signal, that integrates down 
with the number of clusters, and the instrumental noise, that scales with the 
number of pixels; the latter decreases with integration time and is subdomi- 
nant for the Wilkinson Microwave Anisotropy Probe 5-year data. It is proven 
both analytically and numerically that the errors for 5-year WMAP data are 
— 1 5 a/3/ N c \ ustcrs fiK per dipole component. The relevant components of the bulk 
flow velocity are measured with a high statistical significance of up to > 3 — 3.5<r 
for the brighter cluster samples. We discuss different methods to compute error 
bars and demonstrate that they have biases that would over predict the errors, 
as is the case in a recent reanalysis of our earlier results. If the signal is caused 
by systematic effects present in the data, such systematics must have a dipole 
pattern, correlate with cluster X-ray luminosity and be present only at cluster 
positions. Only contributions from the Sunyaev-Zeldovich effect could provide 
such contaminants via several potential effects. We discuss such candidates apart 
from the bulk-motion of the cluster samples and demonstrate that their contri- 
butions to our measurements are negligible. Application of our methods and 
database to the upcoming PLANCK maps, with their large frequency coverage, 
and in particular, the 217GHz channel, will eliminate any such contributions and 
determine better the amplitude, coherence and scale of the flow. 
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1. Introduction. 

Peculiar velocities are deviations from the uniform expansion of the Universe. In the 
gravitational instability model, they are generated by the inhomogeneities in the matter 
distribution. Most determinations of the peculiar velocities are based on surveys of individual 
galaxies. Early measurements by Rubin et al (1976) found peculiar flows of amplitudes ~700 
km/sec. Using the "Fundamental Plane" (FP) relation (Dressier et al 1987, Djorgovski & 
Davis 1987) suggested that elliptical galaxies within ~ 60/i _1 Mpc were streaming at ~ 600 
km/sec with respect to the CMB (Lynden-Bell 1988). Mathewson et al (1992) used the 
Tully-Fisher (1977 - TF) relation and found that this flow of amplitude 600 km/sec does not 
converge until scales much larger than ~ 60/i -1 Mpc, a result that was in agreement with 
a later analysis by Willick (1999). Using the brightest galaxy as a distance indicator for a 
sample of 119 rich clusters, Lauer & Postman (1994 - LP) measured a bulk flow of ~700 
km/sec on a scale of ~150/i _1 Mpc, but a re-analysis of these data by Hudson & Ebeling 
(1997) taking into account the correlation between the luminosity of brightest galaxy and 
that of its host cluster found a reduced bulk flow pointing in a different direction. Using the 
FP relation for early-type galaxies in 56 clusters Hudson et al (1999) found a similar bulk 
flow as LP and on a comparable scale, but in a different direction. A sample of 24 SNIa 
showed no evidence of significant bulk flows out to ~ lOO/i -1 Mpc (Riess et al 1997), and 
a similar conclusion was reached with a TF-based study of spiral galaxies by Courteau et 
al (2000). Kocevski & Ebeling (2006) analyzed the contribution to the peculiar velocity of 
the Local Group due to structures beyond the Great Attractor and found that the dipole 
anisotropy of the all-sky, X-ray-selected cluster sample compiled there suggested that most of 
the flow was due to overdensities at > 150/i~ 1 Mpc. Watkins et al (2009) developed a method 
to suppress the sampling noise in the various galaxy surveys and showed that all the data 
(except for the LP sample) agreed with substantial motion on a scale of ~ 50 — 100/i~ 1 Mpc. 
In a follow-up study Feldman et al (2009) estimated the source of the flow to be at an effective 
distance larger than 200/i _1 Mpc; they suggested that the absence of shear is consistent with 
the attractor being at infinity, as proposed in Kashlinsky et al (2008, hereafter KABKE1). 

Cosmic Microwave Background (CMB) temperature fluctuations in the direction of clus- 
ters of galaxies provide an alternative method to measure peculiar velocities. The scattering 
of the microwave photons by the hot X-ray emitting gas inside clusters induces secondary 
anisotropics (Sunyaev & Zeldovich 1970, 1972) that are redshift independent and, if the 
noise is isolated, can be used to probe the velocity field to much higher redshifts than with 
galaxies. The Sunyaev-Zeldovich (SZ) effect has two components: thermal (tSZ), due to 
thermal motions of electrons in the potential wells of clusters and kinematic (kSZ) due to 
motion of the cluster as a whole with respect to the isotropic CMB rest frame (see review 
by Birkinshaw 1999). However, such measurements for individual clusters are dominated 
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by large errors. On all-sky CMB maps, the bulk flow motion of clusters of galaxies can be 
obtained by using large all-sky cluster samples and evaluating the CMB dipole at cluster 
locations (Kashlinsky & Atrio-Barandela, 2000 - hereafter KAB). We have applied the KAB 
method using the largest - at that time - sample of galaxy clusters in conjunction with the 
3-yr WMAP data and have uncovered a large scale flow of amplitude 600 — lOOOkm/s extend- 
ing to at least ~ 300k- 1 Mpc (Kashlinsky et al 2008, 2009a - KABKE1,2). That analysis 
has now been extended to a still larger and deeper sample of over 1,000 clusters and 5-yr 
WMAP data (Kashlinsky et al 2009b - hereafter KAEEK). The KAEEK analysis confirms 
the KABKE1,2 results and shows that the flow remains coherent and extends to at least 
twice the distance probed in KABKE1,2. A larger cluster sample enabled KAEEK to bin 
the signal by cluster X-ray luminosity (Lx)- The dipoles evaluated for binned subsamples 
increase systematically with increasing L^-threshold, as expected if the signal is produced 
from the kSZ effect by all clusters participating in the same motion, a correlation that would 
not exist if the signal was produced by a rare excursion from noise or primary CMB. 

Upcoming data from both the long integration WMAP data and the PLANCK mission 
will bring more accurate CMB maps. PLANCK data will be particularly important because 
of its wider frequency coverage, finer angular resolution and lower instrument noise. It is 
imperative to identify the prospects and limitations of the applications of the current KAB 
methodology to these future datasets. Also, alternative methods can test the existence of a 
large scale flow (see Itoh, Yahada & Takata, 2009; Zhang, 2010). 

In this paper we present a detailed analysis of the uncertainties affecting the measured 
bulk flow providing the necessary details to support the KAEEK results. In Sec. 2 we briefly 
summarize our method followed by a summary of data and dipole analysis. In Sec. 3 we 
present a theoretical derivation of the error bars, showing when they become dominated by 
cosmic variance of the cosmological CMB residual that remains in the maps after filtering. 
It is demonstrated that for the CMB sky of our Universe and an isotropic all-sky cluster 
catalog the errors in the KABKE/KAEEK dipole measurements are ~ 15^/3/iV cl /iK per 
dipole component. Given the filtering scheme adopted in our studies, these errors cannot 
be reduced much in CMB data with lower instrument noise. Rather the strategy of further 
increasing the signal-to-noise in the measured dipole value must be through increasing the 
number of clusters. Particularly important would be to increase the number of observed 
clusters at the bright end of the cluster luminosity function, where the much larger cluster 
optical depth, r, compensates for the decrease in the abundance of such clusters (KAEEK). In 
Sec 4 we discuss different methods to estimate error bars describing their various biases. Sec. 
5 addresses the overall statistical significance of the measurement. Given that the measured 
dipole increases with the X-ray luminosity threshold, the signal found in KABKE1,2 and 
KAEEK cannot arise in primary CMB, but tSZ dipole contributions can potentially provide 
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confusion to the measurement. Section 6 discusses such possible systematic effects due to 
tSZ and we show there that all are of negligible amplitude and none could have generated 
the measured signal. Finally, in Sec. 7 we present our conclussions. 

2. Methodology and Analysis. 

2.1. KA-B method 

If a cluster at angular position y has the line-of-sight velocity v with respect to the 
CMB, the SZ CMB fluctuation at frequency v at this position will be S u (y) = 5tsz{y)G(v) + 
SkSz(y)H(u), with 5tsz=rTx/T c>ann and 5 kS z=-T(v p /c) cosd, being 6 the angle between the 
cluster peculiar velocity v p and the line of sight. Here G(y) ~ —1.85 to —1.35 and H{y) = 1 
if the thermodynamic CMB temperature is measured over the range of frequencies probed 
by the WMAP data, r is the projected optical depth due to Compton scattering, T x is the 
cluster electron temperature and /cBT , e ,ann=511 KeV. When averaged over many isotropically 
distributed clusters moving with a significant bulk flow with respect to the CMB, the kine- 
matic term generates a dipole contribution that could dominate, enabling a measurement 
the bulk motion Vbuik of the cluster sample. Thus, KAB suggested measuring the dipole 
component of 5 v {y) at N c \ cluster locations on the CMB sky. 

We denote by (do, o,i x , aiy, a>iz) the monopole and three dipole components evaluated over 
some locations in the sky and follow the same conventions as in KABKE1,2 and KAEEK: 
ao = (AT) and an = (ATVij) with % = (x,y,z) and n XtViZ = (sin 9 cos 0, sin 9 sin 0, cos 9) 
are the direction cosines of a vector with angular coordinates (0,(f>). Brackets represent 
averages taken over the cluster population of our catalog. The definition of monopole and 
dipole above follows the convention used in the Healpix remove_dipole routine. The dipole 
power is defined as C\ = Ylm=-i l a im| 2 > where a lm are the three dipole components. With 
our normalization, C^km is such that a coherent motion at velocity 14, u ik would lead to 
Ci,kin = ^cMB( r ) 2 Kuik/ c2 > where T C mb = 2.725K is the present-day CMB temperature. For 
reference, - l((r)/10- 3 )(\/ bulk /100km/sec) piK. 

When the dipole is computed at the position of iV cl clusters, it will have contributions 
from 1) the instrument noise, 2) the tSZ component, 3) the primary cosmological CMB 
fluctuation component from the last-scattering surface, and 4) the various foreground con- 
tributions at WMAP frequencies. The latter can be significant at the K and Ka WMAP 
channels, so we restricted our analyses to the WMAP Channels Q, V and W which have 
negligible foreground contributions. 
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For N c \ 3> 1 the dipole of the observed 5 U becomes 

n. ~ n kSZ + fl tSZ + fl CMB + ffn ° isC (I) 

Prior to any analysis, the CMB dipole due to our motion with respect to the isotropic 
CMB frame is removed from the data. The KSZ effect measures velocities with respect to 
the CMB frame which also is taken to be the frame of the universal expansion. This does 
not change when all-sky dipole or any other £-pole moments are subtracted in the all-sky 
maps. This dipole subtraction removes our peculiar velocity, tWai, contributions down to 
0[(v\ oca _i/c) 2 } contributions to the quadrupole. To check that the latter does not contribute 
to the measurement, we also ran the pipeline subtracting the all-sky quadrupole from the 
original maps and detected only negligible differences in the final results. As shown in KAB, 
in this way the kSZ term can be isolated in eq. 1. 



2.2. Map preparation and analysis 

The process that enabled us to isolate the kSZ term is described in detail in KABKE1,2. 
Briefly: 

(A) An all-sky catalog of X-ray selected galaxy clusters was constructed using available 
X-ray data extending to z ~ 0.3. 

(B) As indicated, we only used WMAP Q, V and W bands, where the foreground 
contamination is smallest. We applied the 3- and 5-yr version of the KpO mask to remove 
those pixels where galactic or point source contributions dominate. Next, to prevent any 
power leakage from the dipole generated by our peculiar velocity, it was removed from 
the pixels that survived the mask. Furthermore, KAEEK explicitly removed dipole and 
quadrupole from the original maps and demonstrated that the quadrupole did not contribute 
to the results. This removes i>i OC ai down to O[(t>i oca i/c) 3 ] contribution to the octupole. 

(C) The cosmological CMB component was removed from the WMAP data using a 
Wiener-type filter, constructed using the ACDM model that best fit the data. It was con- 
structed in order to minimize the difference ((ST — noise)) 2 ). Next, filtered maps were 
constructed using all multipoles with £ > 4 and keeping the same phases as in the original 
maps. Modes with £ < 3 were not included to avoid any possible contributions that could 
be introduced by the alignment of those low order multipoles and also because those modes 
would potentially be the most affected by any hypothetical power leakage. 

(D) The Wiener filter is constructed (and is different) for each DA channel because the 
beam and the noise are different. This prevents inconsistencies and systematic errors that 
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could have been generated if a common filter was applied to the eight channels of different 
noise and resolution. 

(E) In the filtered maps, the monopole and dipole are computed exclusively at the cluster 
positions, using Healpix remove_dipole routine ascribing to each cluster a given circular 
aperture. Due to the variations of the Galactic absorbing column density and ROSAT 
observing strategy, cluster selection function and X-ray properties may vary across the sky 
introducing possible systematics. In KABKEI,2 we used the measured X-ray extent of 
each cluster, Ox and computed the dipole for different apertures, in multiples of Ox and, 
to avoid being dominated by a few very extended nearby clusters like Coma, we introduced 
a cut so the final extent of any cluster was always smaller than 30'. There we computed 
core radii directly from the data and from an L x — r c relation. Analyses using both sets 
gave consistent results, consistent with the X-ray systematic effects not affecting our results 
significantly. More important, variations in the final aperture were already tiny in the 
KABKEI,2 analysis and KAEEK used altogether a fixed aperture were the mean monopole 
vanishes. The KAEEK results are consistent with the previous (KABKE1,2) measurements. 
Fixing the same aperture for all clusters simplifies the statistical analysis and this is the 
approach taken in this article. 

(F) We compute the monopole and dipole for different angular apertures. At small 
apertures (~ 10'), clusters show a clear tSZ decrement, but the amplitude of the signal falls 
off with increasing angular aperture. The final dipole is computed at the aperture where 
the mean monopole of the clusters vanishes. This ensures that the TSZ contribution to the 
measured dipole is negligible and does not confuse the KSZ component. 

(G) Our final result is a dipole measured in units of thermodynamic CMB temperature. 
To translate the three measured dipoles into three velocity components, we need to determine 
the average cluster optical depth to the CMB photons, (r), on the filtered maps. Since 
filtering reduces the intrinsic CMB contribution, it also modifies its optical depth, r. In 
KABKEI,2 we introduced a calibration factor C^ioo that gave the kSZ dipole in /xK of a 
bulk motion of amplitude \4>uik = lOOkm/s. The calibration factor depends both on the filter 
and on the cluster profile. In KABKEI,2 and KAEEK it was estimated using a (5 model and 
the angular X-ray extent of the cluster. 

We defer to Sec. 5 a discussion on the statistical significance of our measurements. We 
emphasize that in the filtered maps we measure monopole and dipole simultaneously. The 
monopole is dominated by the tSZ component and its amplitude sets an upper limit on a*^f 
(see eq. 1), the tSZ dipole due to an inhomogeneous cluster distribution on the sky. We 
found a dipole at cluster positions with a high confidence level and we obtained this dipole 
at the (fixed) cluster aperture when the tSZ monopole component was zero. Since the tSZ 
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component from the clusters vanishes, only a contribution from the kSZ component, due to 
large-scale bulk motion of the cluster sample, remains. 

The main current uncertainty in our method is the calibration, currently parameterized 
with the Ci ; ioo quantity, which generally is a matrix. At present, we do not have enough 
information on the tSZ profile of the clusters in our catalog to increase the accuracy of 
our calibration. Sec. 8 of KABKE2 discusses the issues and points out that we may be 
overestimating the velocity amplitude in the current cluster catalog by ~ 20 — 30%. However, 
this question is significant only insofar as the precise amplitude of the flow velocity in km/sec 
is concerned. PLANCK, with its large frequency coverage, will allow the measurement of 
individual profiles for an important fraction of clusters in our catalog and should enable us 
determine the calibration coefficients more accurately. 



3. Noise and Intrinsic CMB Residual Contributions. 

The KAB method to measure bulk flows using clusters of galaxies as tracers of the 
velocity field requires the intrinsic CMB component to be removed from the data. To this 
end we have designed a filter which minimizes ((ST — noise) 2 ). As we show below this 
filter removes the primary CMB anisotropies down to the fundamental limit imposed by the 
cosmic variance. In Fourier space this filter is expressed as 

\d e \ 2 - CfBj 

where \d(\ 2 = (21 + \ae m \ 2 is the power measured in each Differencing Assembly 

(DA) corrected for the mask sky area, and Cj h B 2 is the power spectrum of the theoretical 
model that best fits the data, convolved with the antenna beam B^ of each DA. Although this 
filter removes much of the intrinsic primary CMB contributions, it leaves a residual CMB 
component since the theoretical model does not reproduce perfectly the data measured at 
our location. This residual will be common to all frequencies and, since it is correlated 
between the various DA's, it limits the accuracy down to which the primary CMB can be 
removed in the KAB method. 

Because of the cosmic variance, the power of the CMB sky at our location Cj/ oc differs 
from the theoretical model Cf 1 and so a residual CMB signal from primary anisotropies is 
left in the filtered maps. To estimate the contribution of noise and the CMB residual to the 
total power on these maps, let ST(h) = Fea>emYem(n) be the temperature anisotropy of 
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the filtered maps expanded in spherical harmonics Ye m . The variance of any filtered map is: 

As indicated, ST(n) contains the cosmological CMB signal and noise, \de\ 2 = Cf oc Bj + Nf. 
The power spectrum at our location differs from the underlying power spectrum by a random 
variable of zero mean and (cosmic) variance A^ = (£ + \)C\ h / f s k yi where f s k y is the fraction 
of the sky covered by the data (Abbot & Wise 1984). Then, due to cosmic variance, Cf oc = 
Cf 1 ± A 1 / 2 . The above limits on Cg bound the range of (jfu, eq. (3), to: 

In this last expression, the variance of the filtered map depends on two components: the 
residual CMB left due to cosmic variance acvju and the noise <tnju, that is not removed by 
the filter. The latter component integrates down with increasing observing time t G b s as t~^[ 2 
and becomes progressively less important in WMAP data with longer integration time. 

We denote by a\ = ±(2q + 1)(A 2 + iV 2 )(Cf + A 9 + iV,)" 1 and let a\t) = £^ =4 a\ be 
the cumulative variance of the residual map. With these definitions, the total variance of the 
filtered map is a'j il = o 2 iX {£ ma ^). For Healpix maps with N si d e = 512 the maximal multipole 
is ^max — 1024 (Gorski et al 2005). In Figure 1 we plot this cumulative contribution of each 
multipole £, <jfu(£), to the total rms of the map. The solid lines represent the mean and 
rms <Jfu(£) of filtered maps of 4,000 realizations of the Ql DA; the shaded area represents 
the dispersion of those realizations, the dot-dashed line is the same quantity but for the 
filtered Ql WMAP 5-year data. The lower dashed lines represent (Jcvjih the residual CMB 
component, and upper dashed line, the total variance of the map [eq (4)]. The dot-dashed 
line also contains any contributions from foreground emissions; the fact that it lies so close 
to the to the region expected from simulating CMB sky implies that foreground emission 
contributions to (jfu are small. Figure 1 clearly shows that for multipoles below £ ~ 200 the 
cumulative variance of the 5- year WMAP maps u 2 {£) is dominated by the residual primary 
CMB signal from the cosmic variance, even though the total variance of the filtered maps is 
dominated by noise. For the Ql WMAP channel, the mean variance of our simulations was 
aj u ~ 2000 (/iK) 2 out of which ~ 200 (/iK) 2 come from the residual primary CMB signal. 

Finally, Figure 1 indicates that our filter removes the intrinsic CMB down to the fun- 
damental limit imposed by cosmic variance. In this sense the filter is close to optimal, since 
it minimizes the errors contributed to our measurements by primary CMB. In principle, 
one can define a more aggressive filter that, together with the intrinsic CMB, also removes 
the noise leaving only the SZ signal. But filtering is not a unitary operation and does not 
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preserve power. Such a filter would then remove an important fraction of the SZ component 
and would probably reduce the overall S/N. In general, a different filter would give different 
dipole (measured in units of temperature) and would require a different calibration. Discus- 
sion of filtering schemes that maximize the S/N ratio and minimize the systematic error on 
the calibration will be given elsewhere. 



4. Monopole and dipole uncertainties. 

Here we consider how the two components present in the filtered maps, i.e. 1) residual 
primary CMB and 2) instrument noise, contribute to the uncertainty in the measurement 
of bulk flows. In KABKEI,2 we adopted two methods to estimate the uncertainties: (I) 
evaluating monopole and dipole on the filtered maps outside cluster locations and (II) using 
the same cluster template on simulated maps. Both methods are different but complemen- 
tary. Errors estimated using method I include any contribution originated by foreground 
residuals and CMB masking while in Method II we account for the inhomogeneity of the 
cluster distribution on the sky. 

It is important to emphasize that the filtered maps have no intrinsic monopole or dipole 
by construction. Since we measure these two moments from a small fraction of the sky, our 
limited sampling generates an error due to (random) distribution of these quantities around 
their mean zero value. The sampling variances of (a ) and (an) are Var((a )) = (afy/N, 
Var((ai)) = (af^/N, where N is the number of independent data points. Direct computation 
shows that: 

a^(al) = ((ATY) a} ^ (a*) = i = (x,y,z) (5) 

In this expression, are the direction cosines of clusters. If clusters were homogeneously 
distributed on the sky then (nf) = 1/3 and one should recover the dipole errors of <7j = y/3a . 
Thus the error on the monopole serves as a consistency check in any such computation. 

Sec. 3 discussed the two components of the variance of the filtered map. As before, 
a cvju an d a N,fu represent the contribution to the total variance due to the residual CMB 
component and the noise, respectively. When we estimate error bars by placing random 
clusters on the real filtered maps outside clusters (Method I) or the real clusters on simulated 
filtered data (Method II), N c i clusters occupy N P i X in No a Differencing Assemblies. As the 
residual CMB signal is correlated from map to map, it will decrease only as the number of 
clusters increases, but the noise term will decrease much faster since it is uncorrelated from 
map to map and pixel to pixel and integrates down with Nda and the integration time. 
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Then, the resulting sampling variance will be 

2 a CVJil , ^NJilitobs) 2 a o 2 

As expected from eq. 5, for an homogeneous cluster catalog the variance in each dipole 
component is three times larger than on the monopole since three quantities are derived 
from the same data set. From Figure 1 we obtain that ocvju — 15//K and (Jnju — 40/iK. 
When clusters are not homogeneously distributed in the sky, the basis of direction cosines is 
no longer orthogonal and error bars need to be estimated numerically. 

The results presented in Figure 1 together with eq. (6) indicate that Method II will 
give slightly larger error bars. If monopole and dipole are evaluated at cluster positions on 
simulated maps then o~cvja an d ^nju in eq. (6) will be close to the average CMB residual 
and noise of the simulated maps. As Figure 1 indicates, they are larger than the filtered 
data (shown by a dot-dashed line) corresponding to the CMB realization representing our 
Universe. The latter, however, is the only CMB sky relevant for the true error analysis in 
this measurement. This fact was already noticed in KABKE2 where such comparison was 
made and the errors were found to be 10-15% larger if using Method II. 

To avoid this bias, we introduce Method Ha: error bars are computed from random 
realizations of the power spectrum of the filtered maps. In Figure 2 we plot the histograms 
of the monopole and dipole components of 4,000 simulations of 1000 clusters with constant 
angular size of 30' with both Method I (random clusters located outside the mask on the real 
data) and Method Ha (the cluster template is fixed and the sky is simulated; the spectrum 
are gaussian realizations of the measured power of the filtered maps). From left to right 
we display the histogram of the monopole and (x, y, z) components of the dipole. The rms 
deviations, given to the left and right of each plot, correspond to Method I and Method Ha, 
respectively. Solid lines represent the histograms in Method I and dashed lines in Method Ha. 
We find that, to good accuracy, the distribution of the monopole and dipoles is Gaussian with 
zero mean. More importantly, we see no systematic differences between both methods. Then, 
neither foreground residuals nor cluster inhomogeneities have a significant contribution to 
the estimated error bars. Instead the errors are dominated by the sampling/ cosmic variance 
when measuring the monopole and dipole from a limited fraction of the sky. 

To test the validity of eq. (6) we carried out another 4,000 simulations with different 
number of clusters: N c \ =100, 180, 320, 570 and 1000 in accordance with Table 1 of KAEEK. 
In Method I, we placed N d clusters at random on the sky. To be fully consistent with how 
cluster samples are selected from the data, the smaller samples are subsets chosen randomly 
from the full sample. In Figure 3 we plot the rms deviation of the monopole (open triangles) 
and the three dipole components. Filled circles, diamonds and solid triangles correspond to 
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the (x, y, z) dipole components. Solid lines connect the results when clusters are assigned 
radii of 30', while dashed lines correspond to results with 20' clusters and follow the same 

1/2 

ordering as the solid lines. The figure shows that 07o,s,y,z) oc NJ with great accuracy. 
As expected, the errors are larger when the cluster size is smaller because of the different 
numbers of pixels entering the instrument noise contribution in eq. (4). In Figure 3, the 
differences between the dipole components come from differences in sky coverage. The x and 
y components, that are in the plane of the Galaxy, are determined with progressively less 
accuracy since the CMB data in the Galactic plane is dominated by foreground emission. 
Still the difference from the uncertainty of the z-component is small (< 10%) particularly 
for the better measured y-component of the dipole. 

We can use eq. (6) to estimate how accurately we measure any dipole component com- 
pared to the monopole. In Figure 4a we plot o^^^ja^ the ratio of the rms deviation of 
dipole components to the rms of the monopole of the 4,000 simulations generated using 
Method I, as described above. Filled circles, diamonds and triangles correspond to the ratio 
of the (x, y, z) rms deviation of the dipole to that of the monopole, respectively. In Figure 4b 
we plot the same magnitudes for Method Ha. The dotted line represents cr( x ,y,z) f^o — V5 
that corresponds (in eq 6) when clusters are homogeneously distributed in the sky. Fig- 
ure 4a clearly indicates that when clusters are chosen randomly on the sky, the error on the 
x-component is larger than on y or z and the scaling with the number of clusters was very 
close to iV*/ 2 , as seen also in Figure 3. In Figure 4b the behavior is very similar: the error 
on the x-component is largest. In this case, and since in Method Ha the cluster template 
is fixed, the scaling is not as exactly oc N^ 2 , reflecting the inhomogeneities present in the 
cluster distribution. However, these deviations are not very significant. 

To study the effect of cluster inhomogeneities potentially present in studies based on 
other catalogs, we carry a different analysis. In Fig 4c we excise clusters from the KAEEK 
catalog as a function of galactic latitude. We plot (n 2 )^ 1 ^ 2 evaluated over the cluster dis- 
tribution when all clusters with \b\ < b cut are removed. Thick and thin solid and dashed 
lines correspond to the x,y,z components, respectively. The dotted line is (h 2 ) _1//2 = \/3, 
that would correspond to a cluster catalog that samples the sky homogeneously. Since the 
mask removes the data in the Galactic Plane, for b cu t < 20° there is little deviation from the 
KAEEK errors. For much larger values of b cut) the error on the x and y components increases 
while the error on the z component approaches that of the monopole. Then, eq (6) permits 
to write the error bars as (J( x ,y,z) — (1-12, 1.05,0.87) x 15-^3/AT cl /iK, i.e., the expected accu- 
racy for each of the components would be only 12 and 5% worse that for an all sky survey, 
compared with that of the monopole, while the z component would be 12% better since the 
Galaxy removes the region of the sky where there is no contribution to it. 
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Comparing the different panels in Figure 4, we see that a y /ao may be smaller than 
the value estimated from the geometry of the catalog as is evident when comparing this 
ratio for N c i = 1000 in Figure 4b with Figure 4c. However, while in (c) the ratio of the 
errors is computed from the cluster geometry, in (b) they are estimated from simulations 
whereby in Method Ha we use simulations of the power spectrum of the filtered CMB data. 
Since monopole and dipole are sensitive to different parity multipoles (even vs odd), the 
slightly lower value of the dipole components with respect to the monopole is reflecting a 
power asymmetry between odd and even multipoles in the filtered map. So, on average the 
monopole is larger than in a random sky and the dipole is smaller. This effect introduces an 
extra variance and enhances the differences between Fig 4a and b. 

When this paper was being completed, Keisler (2009) replicated the analysis of KABKE1,2 
compiling his own X-ray cluster catalog using publicly available data. Analogously to 
KAEEK and this study, he noticed that the errors on WMAP 3-year data were already 
dominated by the residual CMB and not by the noise. He confirmed the measured central 
dipole values of KABKE2, but claimed significantly larger errors than KAEEK, particularly 
for the y-component. (We note again that if the KABKE1,2 dipole originated from primary 
CMB and/or noise, its magnitude should display no correlation with the cluster luminosity 
threshold that was demonstrated to exist in KAEEK). Specifically, in the final configuration 
his catalog contained ~ 700 clusters and his claimed errors were crf^ CT ~ (1.7, 1.7, 1.1) /jK. 
Those errors are larger than those quoted in KAEEK. A small increment (of order of 10- 
15%) can be accounted for by his treatment of the errors using simulations of the CMB sky 
around the theoretical ACDM model and thereby pumping up the cosmic variance compo- 
nent (see Fig. 1), as well as anisotropies in his catalog. Keisler (2009) uses a catalog without 
recomputing cluster properties from X-ray data, a procedure done in Kocevski & Ebeling 
(2006). That dataset is then less complete, especially at low latitudes, but that in itself can 
account only for a small increase in the errors. However, Keisler (2009) claims an increase 
in errors by a factor of > \/20 compared to KABKE2. Clearly, the effect of residual CMB 
correlations between the N D a = 8 WMAP channels can at most increase the KABKE1,2 
errors by a factor of \fN D A < \/8. (In reality, because the instrument noise is also present, 
the errors on individual dipole components in Table 2 of KABKE would be increased for 
3- year WMAP data by a factor of ~ \/G to become < 1/iK at the largest redshift bin.) A 
larger increase, as we have demonstrated above, cannot happen. In our computations we do 
not reproduce Keisler (2009) errors with proper analytical and numerical procedures, even 
using his methodology. 

Interestingly, we recover the magnitude of his claimed errors if one important aspect 
of the KABKE processing is omitted. When working on simulated data such as in Method 
II, care must be taken to replicate all the details of the data analysis done in KABKE1,2. 
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The filter must be constructed using the theoretical model and the simulated data. Since 
only modes with i > 4 are used to generate the filtered map, by construction it has zero 
monopole and dipole. But this map covers the full sky and is not yet the correct model for 
the data. One needs to remove the monopole and dipole outside the Galactic mask, as is 
done with the filtered data in our processing. The full map has zero monopole and dipole, 
but the fraction of the sky outside the mask does not. To test this effect we carry out two 
sets of 4,000 simulations SI and S2, starting with the same initial seed and using Method 
II. In SI simulations we removed monopole and dipole outside the Galactic mask; in S2 we 
did not. In Figure 5 we show the histograms with the distribution (from left to right) of the 
monopole and the (x,y,z) components of the dipole. The solid and dashed line shows the 
results for the SI and S2 simulations, respectively. The labels on the left (right) give the rms 
deviation for SI (S2). The differences can be easily explained: if the monopole and dipole 
are not subtracted the measured monopole and dipole at cluster locations are not different 
from zero simply because we are sampling the signal over a very small fraction of the sky. 
Rather they are not zero because we are measuring the monopole and dipole present on the 
fraction of the sky outside the mask. We checked that when in the SI simulations we add 
the variance of the monopole/dipole subtracted outside the mask and the variance of the 
monopole/dipole computed at cluster positions, we obtain exactly the variance measured 
in the S2 simulations. For instance, the variances on (ao, ai x , ai y , a± z ) outside the mask are 
(0.5, 1.5, 1.3, 0.1)(/iK) 2 . The variances on the monopole/dipoles measured at the location of 
1000 clusters of our catalog are (0.5, 1.8, 1.3, 0.7)(/iK) 2 ; if added with the previous variances, 
the monopole/dipoles error bar increase by (46,36,41,5)%, respectively. As expected we 
see that, since the z axis is perpendicular to the galactic plane, the error bars are boosted 
preferentially in the x and y directions. This explains that in S2 simulations the error in 
the monopole o"o - as well as a X)V - is larger than a z . Only in the case of the faulty S2 
processing do we recover the magnitude of errors found by Keisler. We cannot claim that 
this step was necessarily overlooked by him but we do find this coincidence puzzling especially 
when considering the deviation of his ratios of o\ y jo\ z and <J\ x ja\ y from the (analytically) 
explained ratios (Fig. 4) and Fig. 1. 



5. The Statistical Significance of the "Dark Flow". 

Because of the correlations in the final filtered maps between the eight WMAP DA's, 
the S/N of the Dark Flow measurement is smaller than suggested in Table 2 of KABKE2 
for individual dipole components, although not fatally so. This was corrected in KAEEK, 
where it was also demonstrated that the dipole correlates strongly with the cluster X-ray 
luminosity L x , as it should if the dipole signal originated from the kSZ effect and not from 
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the primary CMB. The minimal S/N of the dark flow measurement is, of course, given by 
the single DA map processing. Fig. 8 of KABKE2, which plots the mean CMB temperature 
decrement over cluster pixels versus the cosine of the angle between the cluster and the apex 
of the motion, shows that KABKE1,2 already detect the dipole at cluster positions at the 
~ (2.5 — 3)<t level in each of the eight DAs. The overall S/N cannot then be lower than this 
floor level. KAEEK further increase this significance and measure the motion to a much 
larger scale. The systematic uncertainties in our calibration procedure do not yet allow us to 
quantify the properties of the flow better, but we hope to accomplish this task in the coming 
years. 

In the KAEEK catalog, the error on the y component is only 5% larger than what 
it would be for a homogeneous cluster catalog. Future versions of the catalog will include 
clusters at higher redshifts that will help to probe the velocity field on even larger scales. 
A great effort is devoted to produce a spatially homogeneous and flux limited sample. If 
the "Dark Flow" is but a large scale flow that affects all the scales out to the horizon, 
one could argue that the signal is uniform on the entire sky and would be unaffected by 
anisotropies on the cluster distribution in alternative catalogs, but this is not so, as Fig. 4 
indicates: incompleteness and asymmetries increase the error bars and could make some 
cluster catalogs insensitive to the flow. 

The original evidence in favor of the measurement being real were three (KABKE1,2): 
(a) the motion was found at cluster positions, (b) it was persistent when the number of 
clusters increased from < 150 to > 700, (c) the dipole kSZ signal was measured when the 
tSZ monopole vanished. Since the thermal and kinematic components are both generated 
by the X-ray gas, it was thought that a measurement of the kSZ effect could be obtained 
only when enough frequency coverage allowed to remove the thermal contribution, because 
of their different frequency dependence. However, in Atrio-Barandela et al (2008 - hereafter 
AKKE) we showed - for the first time - that cluster gas distribution follows an NFW profile 
(Navarro, Frenk & White 1996). Then, cluster temperature falls with radius and, by adding 
the contribution from the cluster outskirts, the kinematic component dominates over the 
thermal in the KAB method (KABKE2). If clusters were isothermal, the thermal SZ signal 
dipole due to the inhomogeneous distribution of the sky could be large enough to make the 
kSZ effect undetectable at WMAP frequencies. 

In KAEEK we provided further evidence in support of the cluster bulk flow being a real 
effect. The cluster catalog used there was large enough to allow the analysis to be carried 
out in luminosity bins. The kSZ signal is AT ksz ~ tvb- Since r is proportional to the cluster 
electron density, it correlates with X-ray luminosity. If the velocity does not correlate with 
cluster luminosity, for example if the flow is homogeneous across the cluster sample, we 
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would expect the dipole evaluated at different cluster subsamples to be larger for the more 
luminous clusters. In KAEEK we were able to carry out such test by decomposing the 
sample in luminosity bins and the analysis conclusively showed that (d) the measured dipole 
correlates with X-ray luminosity, strengthening the evidence against a possible undiagnosed 
systematic effect. 

In KAEEK it was shown that clusters with the highest luminosity dominate the S/N 
of the measured flow. To quantify the level of statistical significance there, we generate 
10,000 dipole components drawn from a gaussian distribution with zero mean and rms the 
measured error of each component as shown in Table 1 of KAEEK. The significance is then 
the percentage of simulated values that deviate from zero less than the measurement. For 
instance, when we consider the measured dipoles for L x > 2 x 10 44 erg/sec clusters with 
z < 0.25 we measured (a lx: a ly , a lz ) = (3.7 ± 1.8,-4.1 ± 1.5,4.1 ± 1.5)/xK. If the dipoles 
an are Gaussian-distributed random variables, the amplitude of the flow for these clusters 
is detected at the 99.95% level consistent with our simulations (in Method I we find just 2 
realizations out of 4,000 with such parameters). For some other configurations in Table 1 
of KAEEK the confidence level would be even higher. Foreground contributions, by their 
non-Gaussian nature, can in principle alter the above percentiles, but the fact that our 
Universe lies so close to the lines in Fig. 1 generated from pure primary CMB, implies 
that foreground emissions contributions are small in our calculations. We do not necessarily 
advocate the above levels to be highly precise, but this discussion clearly shows that we 
recover a very statistically significant dipole. While the dipole components are less significant 
in lower L^-bins, presumably because of the lower r's for these clusters, the a ly component is 
always negative and a\ z almost always positive in all three L^-bins, while the a\ x component 
oscillates and is the least accurately measured component. In this case, the possibility that 
a± x is zero can be rejected at more than 95% and a\ z ,a\ z at the 99% confidence level. Due 
to the changing sign, the measurements of the lower Lx bins reduces the significance of the 
detection of a lx and we can not claim any measurement but the other two components are 
still significant at more than 95%. Finally, these probabilities would become even higher if 
one folds in the directional coincidence of the recovered dipole to that measured by Watkins 
et al (2009) from galaxy surveys data on smaller scales, < 100 Mpc. 



6. Possible, but negligible, Lx-dependent (SZ) systematics 

The only possible systematic effect that could mimic our measurements would have to 
be present exclusively at cluster positions, produce zero monopole and also give a dipole 
which increases with increasing Lx- Such systematics cannot come from primary CMB, and 
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would have to originate from contributions by the SZ components, which depend on Lx 
in the appropriate manner. Since in KAEEK we showed that the measured dipole corre- 
lates with the X-ray luminosity threshold, it is important to discuss possible L^-dependent 
contributions even if only to rule them out because of their negligible magnitudes. Given 
that we evaluate the dipole at the aperture where the monopole vanishes, there are three 
ways that could potentially confuse the measurement: 1) Systematic effects that could fold 
the Doppler-shifting due to the local motion into the tSZ contributions, 2) cross-talk effects 
between the tSZ monopole and dipole terms in sparse/small samples (Watkins & Feldman 
1995); and 3) inner motions of the intracluster medium (ICM) as opposed to the coherent 
flow of the entire cluster sample. 

We discuss all three of these contributions below and demonstrate that they are negligi- 
ble. Before we go into the rest of the section, we emphasize again that the dipole at cluster 
positions is measured at zero monopole. That monopole vanishes within the noise with 1-a 
uncertainty of ~ 15/v^/Vd/iK or amplitudes significantly below 1/iK for N c \ > 200; the actual 
numbers are given in Table 1 of KAEEK. 

The first two of these contributions come from the tSZ component, while the latter would 
arise from the kSZ effect. Thermal and kinematic SZ dipoles will differ in one very important 
aspect: their frequency dependence. WMAP measures only in the Rayleigh- Jeans part of 
the spectrum and for the Q, V and W bands the change in amplitude is about < 30%, too 
small to be distinguished (see KABKE2 for a discussion). A tSZ induced dipole will change 
sign in the Wien part of the spectrum, while a kSZ dipole will preserve it. The latter will 
be different from zero at 217GHz, the zero crossing frequency of the thermal component. 
Although we show below that the tSZ induced contributions to the dipole are very small, 
PLANCK with its large frequency coverage covering both sides of 217 GHz will be definitive 
in this respect. 



6.1. Systematics due to tSZ shift from the local motion 

The intrinsic CMB dipole due to the motion of the Sun is over two orders of magnitude 
larger than the measured cluster dipole. This motion is known to be u & — 370 ± 3km/s in 
the direction (l,b) = (264°, 48°), close to the direction (276°, 30°) of the Local Group with 
respect to the same reference frame (Kogut et al, 1993) and is not far within the errors 
from the direction measured in KAEEK: (290 ± 20, 30 ± 15)°. An undiagnosed systematic 
effect, present in the time ordered data or in our pipeline that affect preferentially the tSZ 
signal, could fold the motion of the Sun into our measurement. For example, a residual 
of the CMB all-sky dipole (AT) res coupled to the thermal SZ effect would correlate with 
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X-ray luminosity and would satisfy the same properties (a-d) as the kSZ effect, except its 
frequency dependence. The amplitude of such undiagnosed systematic dipole will be bound 
by (AT) res < (AT) t sz(uo/c). In AKKE we showed that the tSZ amplitude of clusters in 
unfiltered maps is of the order of ~ — 30/iK and this amplitude is reduced a factor of ~ 3 
due to filtering (KABKE2). Then, any possible systematic effect that correlates with cluster 
luminosity will be (AT) res < 10 -2 /iK, more than 2 orders of magnitude smaller than the 
measured effect. 



6.2. Cross-talk from tSZ monopole in KAEEK sample 

Since clusters are not randomly distributed on the sky, the tSZ signal will give rise 
to a non-trivial dipole signature that, in principle, may confuse the kSZ dipole. The tSZ 
dipole for a random cluster distribution is given by ~ ((AT) tS z)(3/AT c i)- 1 / 2 decreasing 
with increasing N d . This decrease could be altered if clusters are not distributed randomly 
and there may be some cross-talk between the monopole and dipole terms especially for 
small/sparse samples (Watkins & Feldman 1995). As discussed in KABKE2, the dipole 
from the tSZ component varies with the cluster sub-sample, contrary to measurements, and 
also has negligible amplitude because it is bound from above by the remaining monopole 
amplitude of ((AT) tSZ ) <C 1/iK measured at the final aperture (see Table 1 of KAEEK). 

In order to assess that there is no cross-talk between the remaining monopole and dipole 
which may confuse the measured kSZ dipole, we proceed in the same manner as in KABKE2 
(see Fig 6 there) repeating the following experiment: 1) The tSZ and kSZ components from 
the catalog clusters were modeled using cluster parameters derived for our current catalog. 
To exaggerate the effect of the cross-talk from the tSZ component, the latter was normalized 
to ((AT)tsz) = —1/iK, a value significantly larger than the monopoles in Table 1 of KAEEK 
at which the final dipole was measured; the results for even larger monopoles were also 
computed and can be scaled as described below. For the kSZ component each cluster was 
given a bulk velocity, Vbuik, in the direction specified in Table 1 of KAEEK, whose amplitude 
varied from to 2,000 km/sec in 21 increments of 100 km/sec. The resultant CMB map 
was then filtered and the CMB dipole, ai m (cat), over the cluster pixels computed for each 
value of Vbuik- 2) At the second stage we randomized cluster positions with (I, b) uniformly 
distributed on celestial sphere over the full sky for a net of 500 realizations for each value of 
Vbuik- This random catalog keeps the same cluster parameters, but the cluster distribution 
now occupies the full sky (there is now no mask) and on average does not have the same 
levels of anisotropy as the original catalog. We then assigned each cluster the same bulk flow 
and computed the resultant CMB dipole, ai m (sim), for each realization. The final ai m (sim) 
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were averaged and their standard deviation evaluated. 

Fig. 6 shows the comparison between the two dipoles for each value of \4uik for the 
most sparse sub-samples from Table 1 of KAEEK. We also made the computations at tSZ 
monopole values still larger than above (see upper left panel for one such example). The 
overall contribution from the tSZ component to the dipole is oc (AT t sz), so in the absence 
of cross-talk effects the amplitude of the scatter in the simulated dipoles is made of two 
components: 1) remaining tSZ oc (AT t sz) and 2) genuine kSZ dipole with amplitude oc Vbuik 
to within the calibration. One can see that there is no significant offset in the CMB dipole 
produced by either the mask or the cluster true sky distribution. The two sets of dipole 
coefficients are both linearly proportional to Vbuik and to each other; in the absence of any 
bulk motion we recover to a good accuracy the small value of the tSZ dipole marked with 
filled circles. As discussed in KABKE2, since the bulk flow motion is fixed in direction and 
the cluster distribution is random, one expects the calibration parameterized by C^ioo to be 
different from one realization to the next, e.g. in some realizations certain clusters may be 
more heavily concentrated in a plane perpendicular to the bulk flow motion and the measured 
Ci.ioo would be smaller. In our case, the mean C^ioo differs by < 10% suggesting that our 
catalog cluster distribution is close to the mean cluster distribution in the simulations. This 
difference in the overall normalization would only affect our translation of the dipole in /zK 
into Vbuik in km/sec, but we note again the systematic bias in the calibration resulting from 
our current catalog modeling clusters as isothermal /3-model systems rather than the NFW 
profiles required by our observations (AKKE, KABKE2). We have no progress to report 
on this issue beyond discussion in sec. 8 of KABKE2 and this paper does not address the 
measured velocity amplitude stemming from calibration; this work is in progress and will be 
addressed after the recalibration of our catalog has been successfully completed. 



6.3. Contribution from intracluster flows 

The intra-cluster medium (ICM) may not be at rest in the cluster potential wells as 
a result of mergers during cluster formation process. In principle, our measurement and 
interpretation then may be affected by turbulent motions that give rise to a kSZ effect that 
would be larger for the more massive clusters. However, since the motions are randomly 
oriented with respect to the line of sight, they will not produce a significant effect. In order, 
to reach the value comparable to \4 u ik ~ 1,000 km/sec, a typical cluster in our sub-sample 

1/2 

of N c \ would need to have thermal motions of ~ Vb u \kN c { , over an order of magnitude larger 
than the velocity dispersion of Coma-type clusters. Rather these motions will enter the 
overall dispersion budget (noise, gravitational instability and this component) around the 
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coherent bulk flow component shown in Fig. 2 of KAEEK. 

7. Conclusions. 

We have analyzed the statistical significance of the results presented in KAEEK. We 
have identified the main contributions to the error budget: noise and the residual CMB 
contribution. While the instrument noise was important in WMAP 1-year data, it was 
much less so in the 3- and 5-year data. With our filtering scheme, there remains a residual 
contribution due to cosmic variance, which correlates at different frequencies and decreases 
only as the number of clusters increases. We have discussed methods to compute the errors 
and presented analytical discussion to estimate the various contributions to the final error 
budget. Measuring dipoles with a fixed template over simulated skies increases the error bar 
in two respects: clusters do not sample the sky homogeneously and different maps will have 
different CMB residuals. Since the measured CMB sky in our Universe has less power than 
the average ACDM realization, this also can boost the errors, but by only ~10-15%. Also, the 
inhomogeneities on the cluster distribution make the error on the various dipole components 
different. For the y component the increment is about 5% compared to the ideal case. 
We have argued that a proper method to compute error bars would be to perform random 
simulations of the measured power of the filtered maps corresponding to the CMB sky of 
our Universe. However, we found that the difference with taking random clusters outside 
the mask, but using real data, was insignificant. 

We have discussed the evidence supporting the existence of the Dark Flow. Indepen- 
dently, different groups using galaxies as tracers of the density or velocity field are showing 
the amplitude and direction of the local flow that are consistent, albeit at a much smaller 
scale, with the Dark Flow motion (Kocevski & Ebeling 2006, Watkins et al 2009, Feldman 
et al and 2009, Lavaux et al 2009). This analysis with the forthcoming PLANCK data will 
provide an important consistency check. With a scanning strategy different from WMAP 
and with better frequency coverage, it will permit us to characterize still better any possible 
undiagnosed systematic. The 217 GHz band with ~ 5' resolution will be specially useful 
since it will allow to measure the kSZ signal from central parts of the clusters in our catalog 
uncontaminated by the thermal component. 
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Fig. 1. — Cumulative rms deviation as a function of multipole. Solid line and shaded area show 
mean and rms of 4,000 simulated Ql filtered maps. Dashed lines represent the residual CMB 
component of the filtered maps due to cosmic variance, computed using eq. (4) and the residual 
CMB plus the noise components. The dot-dashed line corresponds to the actual Ql band of WMAP 
5-yr data. 
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Fig. 2. — Histograms of the distribution of monopoles and the three dipole components computed 
using the filtered Ql WMAP 5-yr map data. Solid, dashed lines correspond to Method I and 
Method Ha of 4,000 simulations (see text), respectively. Also indicated is the rms dispersion (in 
micro Kelvin) for Method I (left) and Method Ha (right). 
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Fig. 3. — (a) Rms deviation of the monopole and three dipole components computed scaled by the 
number of clusters. Open triangles, circles, diamonds (blue) and solid triangles (red) correspond to 
the monopole and (x,y,z) components of the dipole. Solid lines joint the symbols of clusters with 
30' radius, while dashed lines follow the same ordering than solid lines but correspond to clusters 
with 20' radius. 
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Fig. 4. — Dipole to monopole error bar ratio, (a) (Black) circles, (blue) diamonds and (red) 
circles correspond the the ratio of the (x,y,z) component of the dipole to the monopole, respectively. 
Monopole and dipole were computed using Method I. (b) Same as (a) but monopole and dipoles 
are computed using Method Ha. (c) Ratio of the dipoles to monopole error bars for our cluster 
catalog. The horizontal axis, b cu t indicates that clusters with \b\ < b cu t are excised from the catalog. 
In all three plots, the dotted line represents the ratio for a perfectly isotropic cluster catalog. 



25 




-4 -2 2 4 -2 2 4 -2 2 4 -2 2 4 

T///K T/VK T/ytiK T/jtiK 



Fig. 5. — Histograms of 4,000 realizations of the CMB sky using Ql DA parameters. Solid, dashed 
lines correspond SI, S2 simulations; in SI (S2) the monopole and dipole outside the mask are (are 
not) removed. The left, right rms dispersion corresponds to SI, S2, respectively. The figures are 
given in micro Kelvin. 
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Fig. 6. — The dipole coefficients for simulated cluster distribution (random and, on average, 
isotropic) are compared to that from the true catalog. (See text for details). Each cluster in each 
catalog is given bulk flow of Vbulk from to 2,000 km/sec in increments of 100 km/sec towards the 
apex of the motion from Table 1 of KAEEK. The results from 500 simulated catalog realizations 
were averaged and their standard deviation is shown in the vertical axis. Dotted lines mark the zero 
dipole axis of the panels. The four most sparse samples from Table 1 of KAEEK are shown which 
correspond to the largest Lx-bins giving the best measured S/N. Black/blue/red colors show the 
x/y/z components of a\ m . Filled circles of the corresponding colors show the dipole components 
due to the modeled tSZ component. The upper left panel shows the results for two values of 
the monopoles: in the case of (ATtsz) = — 3^K_,the results are shown as individual error bars; 
the case of (ATtsz) = — is shown with filled contours. All other panels show the results for 
(ATtsz) = — 1//K and our simulations find good scaling with higher monopole values as described 
in the text. 



